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Abstract 

We investigate a special case of the maximum quadratic assignment problem 
where one matrix is a product matrix and the other matrix is the distance matrix 
of a one-dimensional point set. We show that this special case, which we call the 
Wiener maximum quadratic assignment problem, is NP-hard in the ordinary sense 
and solvable in pseudo-polynomial time. 

Our approach also yields a polynomial time solution for the following problem 
from chemical graph theory: Find a tree that maximizes the Wiener index among 
all trees with a prescribed degree sequence. This settles an open problem from the 
literature. 

Keywords. Combinatorial optimization; computational complexity; graph theory; 
degree sequence; Wiener index. 



1 Introduction 

The Quadratic Assignment Problem (QAP) in Koopmans-Beckmann form [1] takes as 
input two n X n square matrices A = (oij) and B = (bij) with real entries, and asks to 
find a permutation vr that minimizes (or maximizes) the objective function 

n n 

Z-n{A,B) := a.T,(^i-)T,^j)bij. (1) 

i=l j=l 
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Here vr ranges over the set 5„ of all permutations of {1,2,..., n}. The QAP is a hard 
and well-studied problem in combinatorial optimization; we refer the reader to the 
book [2] by Qela and the recent book by Burkard et al. [3j for more information on this 
problem. One branch of research on the QAP concentrates on the algorithmic behavior 
of strongly-structured special cases; see for instance Burkard & al [3] or Deineko & 
Woeginger [5j for typical results in this direction. In this paper, we will contribute a 
new topic to this research branch. 

1.1 The Wiener Max-QAP 

This special case of the QAP restricts matrix ^ to be a symmetric product matrix, which 
means that there are non-negative integers ai < a2 < • • • < On such that 

o-ij = otiaj for i < i, j < n. (2) 

The second matrix in the Wiener Max-QAP is the distance matrix of a one-dimensional 
point set, which means that there are integers /3i < /32 < • • • < /3n such that 

bij = for 1 <i,j <n. (3) 

Throughout this paper, a matrix B of the form ([3]) will be called a ID-distance matrix. 
The goal in the Wiener Max-QAP is to maximize the objective value in ([1]), which can 
be rewritten as 

n n 

Z^{A,B) = a^(i)a^(j) Ift (4) 

i=i j=i 

In this paper, we fully determine the computational complexity of the Wiener Max- 
QAP. Our results are as follows. As a negative result, we will prove in Section [2] that the 
Wiener Max-QAP is NP-hard in the ordinary sense. On the positive side, in Section [3] 
we will derive a useful decomposition property, and we will show that there always exists 
an optimal permutation that is V-shaped. These positive results are then applied in 
Section m to derive a pseudo-polynomial time algorithm for the Wiener Max-QAP. 

1.2 The Wiener index of a tree 

The Wiener index W{G) of a connected undirected graph G = {V,E) is the sum of the 
distances between all pairs of vertices in V. The Wiener index was introduced in 1947 
by Harold Wiener [6] to characterize certain molecular structure properties of saturated 
hydrocarbons. We refer the reader to the survey article [7] of Dobrynin, Entringer & 
Gutman for comprehensive information on this fundamental graph parameter. Chemists 
are often interested in the Wiener index of certain trees, where the vertices represent 
atoms and where the vertex degrees correspond to the valencies of the atoms. Entringer, 
Jackson &: Snyder [8] show that among all r-vertex trees, the path Pj. has the largest and 
the star Ki^r has the smallest Wiener index. Fischermann &: al [9] characterize the trees 
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that minimize the Wiener index among all trees with r vertices and maximum degree 
at most A, and they also provide several results on the corresponding maximization 
question. 

Wang [To] describes a simple greedy algorithm for finding a tree that minimizes 
the Wiener index among all trees with a prescribed degree sequence. Zhang &; al [11] 
derive the same result independently and by different techniques. The corresponding 
maximization question remains open (note that the paper [LQl claims a solution to it, 
and that the corrigendum [12] points out a crucial mistake that invalidates these claims) . 
Wang [12] writes about the maximization question: "While the extremal trees seem to 
be difficult to find and are not unique, an algorithm to find at least one of such trees 
may exist and may be easier to find. " 

In this paper, we resolve the computational complexity of the maximization ques- 
tion: there exists a polynomial time algorithm that finds a tree that maximizes the 
Wiener index among all trees with a prescribed degree sequence. More precisely, we 
will show in Section [5] that this maximization question can be modeled as a special 
case of the Wiener Max-QAP (modulo certain minor modifications). Consequently, the 
machinery developed in Sections [3] and S] can be applied to it. The pseudo-polynomial 
time complexity of our algorithm turns into a polynomial time complexity, since all the 
involved numbers are moderately small. 

2 Complexity of the Wiener Max-QAP 

In this section, we establish NP-hardness of the Wiener Max-QAP. The proof is done 
by means of a reduction from the following variant of the partition problem (see Garey 
& Johnson [13]) which is well-known to be NP-hard in the ordinary sense. 

Problem: Partition 

Input: A sequence qi, . . . , q2k of 2k positive integers with Ylfti 1i — '^Q- 
Question: Does there exist / C {1, . . . , 2k} with |/| = k and ^j^j qi = Q? 

We construct an instance of the Wiener Max-QAP of dimension n = 2k. The n x n 
product matrix A is defined by = qi, \/i = 1,2, ... ,n, and hence satisfies Oij = qiqj. 
The ID-distance matrix B uses the points = 1 for \ < i < k and /3j = 2 for 
A: -|- 1 < i < 2k. Note that bij = — /3j| equals 1 if one of the indices i,j lies in the 
range 1, . . . , fc whereas the other index lies in -|- 1, . . . , 2A;; in all other cases bij = 0. 

Consider a permutation tt G S2fc, and let J denote the set of all i with ■n{i) < k. 
Note that |J| = k. By setting x = J2ieJli^ objective value in ([4]) can be rewritten 
as 

XIX] = X^^X^J = x{2Q-x) = 2Qx-x^. 
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As the concave function /(x) = 2Qx — is maximized at x = Q, it is easily seen that 
the QAP has objective value at least if and only if the Partition instance has a 
positive answer. This yields the following result. 

Theorem 2.1 The Wiener Max-QAP is NP-hard in the ordinary sense. □ 



We first discuss a useful decomposition property of the Wiener Max-QAP. Consider 
some fixed instance with product matrix A and ID-distance matrix B, and let I C 
{!,..., n} with |/| = A; be some fixed subset of the indices. Instead of optimizing 
over all possible permutations in ([4]), we only allow permutations vr that map / into 
J = {1, . . . ,k} and that consequently map the complement of / into the complement 
+ 1, . . . , n} of J. Intuitively speaking, this subdivided version QAP assigns the values 
with i G I to the points /3i, . . . , /3k and the values ai with i ^ I to the remaining 
points Pk+i, . . . ,/?„. 

Any permutation vr € 5'„ induces a bijection a that maps J into /, and another 
bijection r that maps the complement of J into the complement of /. Furthermore, 
denote X = Yliei'^i ^ = J2i<^i'^i- The objective value Zt^{A, B) in ([4j) can then 
be written as Zi + Z2 + + Z4 in the following way. 



3 Structure of the Wiener Max-QAP 





(6) 




ieJ j<^.j 



ieJ j^j 




(7) 




(8) 



4 



Note that in the resulting summations in ©-([T]), every single term does either 
depend on function a or on function r, but does never depend on both functions simul- 
taneously. If we collect all the terms in Zi + Z2 + Z3 + Z4 that solely depend on this 
function a, we get 

"-rWaaO') Ift - /3j| + 2^ a„(^i)Y{(3k - ^i). (9) 

We observe that the objective function in ([9]) essentially corresponds to a smaller + 
dimensional Wiener Max-QAP: the underlying ID-distance matrix is built around the k 
points Pi, . . . , /3k plus a duplicated point at (3^ (which corresponds to the occurrence of 
f3k in the right hand sum) . The underlying product matrix is built around the k numbers 
ttj with i G I plus the number Y (which corresponds to the factor Y in the right hand 
sum). Furthermore, the right hand sum in ([9]) imposes the additional restriction that 
the new value Y has to be assigned to the duplicated point at 

As a consequence of all this, the problem of finding the optimal function a and the 
problem of finding the optimal function r are two separate optimization problems that 
can be solved independently of each other. We call this the decomposition property of 
the Wiener Max-QAP. This decomposition property plays a central role in many of our 
arguments. As a first application of the decomposition property, we next deduce a result 
on the combinatorial structure of optimal permutations for the Wiener Max-QAP. 

Definition 3.1 A permutation tt G Sn is called V-shaped, if there exists an index I 
with 1 < i < n such that 7r(i) > 7r(« + 1) for i = 1, . . . and such that TT(i) < 7r(z + 1) 

for i = £,... ,n — 1. 

In other words, a V-shaped permutation tt is first decreasing up to i, and then increasing 
from i onwards, where the increasing or the decreasing part can also be empty. 

Theorem 3.2 Every instance of the Wiener Max-QAP possesses an optimal solution 
TT that is V-shaped. 

Proof. For simplicity of presentation, we will assume without much loss of generality 
that all a-values are pairwise distinct, and that also all /3- values are pairwise distinct. 
The statement for the general case then follows easily from this (by locally reordering 
or renaming the values). 

Now consider an optimal permutation vr, and note that the values vr(l), . . . , 7r(n) and 
a^(i), . . . , a7r(„) are ordered in the same way. Suppose for the sake of contradiction that 
permutation tt has a local maximum at k with ce.^(^k-i) < '^7r(fc) ^^d 0,^(^+1) < Oi^(k)- 
By the decomposition property of the Wiener Max-QAP, the optimal permutation vr 
induces an optimal solution to the five-dimensional problem of assigning the five values 

k-2 n 

L = '^a^{i), a7r(fc-l)i "7r(fe); «7r(fc+l)5 = 

i=l i=k+2 
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to the five points /Sfc-i, l^k-i, l^k, f^k+i, and subject to the constraint that value L 
is assigned to point and that value R is assigned to point /3k+i- 

Now let us switch the positions of a7r(fc-i) Oin{k) ™ the solution that vr induces for 
the five-dimensional problem, such that a7r(fc-i) goes to /3fc and a^^t^k-j goes to (3k-i- Since 
this switch cannot increase the objective value, the difference between the corresponding 
two objective values is non-positive: 

{an(k) - an{k-i)) {l3k - Pk~i) {an{k+i) + R - L) <0 (10) 

In an analogous way, we can switch the positions of a7r(A;) and a^(k+i) iii the induced 
solution. This then leads to the following inequality: 

(a^(fc) - a^(fc+i)) (/3fc+i - I3k) (a7r(fc-i) - R + L) < (11) 

The first two factors on the left-hand side of (jlOh and also the first two factors on 
the left-hand side of (jlip all are positive. This implies a7r(fc+i) + R — L < and 
ct7r(fc-i) — R + L < 0. By summing these two inequalities we arrive at the contradiction 
Q^7r(fc+i) + Q^7r(fc-i) ^ 0. We conclude that vr cannot have any local maximum, and this 
yields that vr indeed is V-shaped. □ 



4 An algorithm for the Wiener Max-QAP 

In this section, we design a pseudo-polynomial time algorithm for the Wiener Max-QAP 
that is based on a standard dynamic programming approach. We recall that the two 
underlying sequences ai < 02 < • • • < On and /3i < /32 < • • • < /?« are in non-decreasing 
order. 

Every state in the dynamic program is specified by a quadruple {k, m, L, R) of inte- 
gers that satisfy the following conditions: 

n 

l<k<n, l<m<n-k + l, 0<L,R, L + R = ^ Ui. (12) 

i=k+l 

With every such state {k, m, L, R) we associate the following {k + 2)-dimensional 
Wiener Max-QAP: the product matrix results from the k + 2 non- negative integers in 
the sequence ai, . . . ,ak, L, R, and the ID-distance matrix results from the k points 
Pm, Pm+i, ■ ■ ■ , Pm+k-1 plus another point in /3m plus another point in /3m+k-i- The goal 
is to find the best solution to this QAP subject to the constraint that the value L is 
assigned to point (3rn and that the value R is assigned to /J^+fc-i- In other words, we 
want to find a bijection a from {m, . . . ,m + k — 1} to {1, . . . ,k} that maximizes the 
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objective value 

m+fc— 1 m+k—l 

^ = X] X] (^<T{i)aa{j)\^i - + 2LR - Pm\ + 

i=m j='m 

m+fc— 1 m+fc— 1 

+ 2 ^ a^(i)L|/3i -/3m| + 2 ^ a„(^i^R\p.m+k-i - Pi\ 

i=m i=m 

We use m, L, i?) to denote the maximum objective value of the corresponding state. 
Next, we will describe how to compute and to store all the values Z{k, m, L, R) step by 
step and in increasing order of k. For k = 1 the corresponding instances are trivial to 
solve, since they only have a single feasible solution. 

Next consider some fixed state {k,m,L,R) with k > 2, and denote M = X^t/*^*- 
By the decomposition property discussed in Section [3] and by Theorem 13.21 there exists 
an optimal bijection a that induces a V-shaped assignment of the k values ai, . . . ,ak 
to the k points /3m, /3m+ij • • • , /3m+fc-i- This implies that Ok as the largest of the k 
values must be assigned either to point /3m or to point /3m+fc-i- First consider the 
case where is assigned to point /?„• Then the remaining k — 1 values are assigned 
to Pm+i, • • • J Pm+k-i, and by the decomposition property the largest possible objective 
value in this case is 

Zi := Z{k-l,m + l,L + ak,R) + 2{L + at) {M + R) \(3m+i - (13) 

In the second case the value Ok is assigned to point fim+k-i- Then the remaining k — 1 
values are assigned to /3m, • • • ,/3m+k-2, and the largest possible objective value is 

Z2 := Z{k-l,m,L,R + ak) + 2{L + M) {R + ak) \/3m+k-i - Pm+k-2\. (14) 

This yields Z(k,m, L, R) = max{Zi, Z2}, and in this fashion one easily determines all 
the function values Z(k, m, L, R) with 2 < k < n. In the end, the optimal objective 
value of the underlying QAP instance can be found as Z(n, 1,0,0). 

Theorem 4.1 The Wiener Max-QAP has a pseudo-polynomial time solution algorithm 
with time complexity 0{n'^ • ^^ai). 

Proof. The correctness of the dynamic programming approach is clear from the above 
considerations. It remains to analyze the time complexity. We observe that there are 
only 0(n^ • X o^i) different states {k,m, L, R) in the dynamic program: there are 0{n) 
possible values for k and m, respectively, and there are 0{Y^ai) possible values for L; 
note that the value of R is already fully determined by the values of k and L. 

The expressions in (jlSp and (jl4p can be evaluated in constant time 0{1): they refer 
to values Z{k — 1, *, *, *) that are known from earlier stages of the dynamic program, 
and they refer to the values M = J2i=i ^^^^ can all be precomputed and stored in a 
preprocessing phase. All in all, this yields that the time complexity is proportional to 
the number of states and hence is 0(n^ • X Q^i)- ^ 
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Note that our approach only yields the optimal objective value. By storing appropri- 
ate auxiliary information in the states of the dynamic program, one can also compute 
the corresponding optimal permutation within the same time complexity. These are 
standard techniques, and we do not elaborate on them. 

5 Maximizing the Wiener index of a tree 

We now return to the Wiener index of a graph that has been introduced and discussed 
in Section 11.21 We will investigate the following algorithmic problem MaxWiener- 
Tree that was left open by Wang \10\ [T2] : an instance consists of a degree sequence 
di, . . . ,dr of r positive integers with Yll=i di = "^^r — 2. The goal is to determine the 
largest possible Wiener index over all trees with degree sequence di, . . . ,dr. 

Recall that a caterpillar is a tree that turns into a path (the so-called backbone of 
the caterpillar) if all its leaves are removed. Shi |14) proved that for every instance of 
MaxWiener-Tree, all the maximizing trees are caterpillars; other proofs of this result 
can be found for instance in Wang [10] and Schmuck [15] . 

Now consider such a caterpillar T, let vi,. . . ,Vn denote the vertices ordered along 
the backbone of the caterpillar, and let £i with 1 < i < n denote the number of leaves 
adjacent to vertex vi. We define set Ci to consist of vertex Vi together with its ii adjacent 
leaves. Then the vertex pairs inside Ci contribute an amount of if to the Wiener index 
W{T), and the pairs with one vertex in Ci and one vertex in Cj [i ^ j) contribute 

£i£j{\j-i\ + 2) + £i{\j-i\ + l) + £^.(|j_i| + i) + 

= {£i + l)i£j + l)\j -i\ + {2£i£j+£i + ej). 

Hence the Wiener index of this caterpillar T is 

n n n 

W{T) = + E E [{£i + '^){£j + l)\j-i\ + {2£i£, +£, + £,)] 

i=l 1=1 j=i+l 

(n \ 2 n ^ n n 

Y,£^] +(n-l)5]£. + -5^5](£. + l)(^, + l)|j-.l. (15) 
1=1 / 1=1 i=l j=l 

Neither the first nor the second sum in ()15p do depend on the way how the numbers 
£i, . . . ,£n are assigned to the backbone vertices vi, . . . ,Vn- The assignment in the third 
sum in (|15p yields an instance of the Wiener Max-QAP with values = £i + 1 and 
points I3i = i. 

Observation 5.1 The problem of finding the maximizing caterpillar for an instance 
o/ MaxWiener-Tree with an explicitly specified backbone vi,. . . ,Vn and an explicitly 
specified sequence £i, . . . ,in of leaf-numbers is equivalent to a Wiener Max-QAP. □ 
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Next, consider a degree sequence di, . . . ,dr that forms an instance of MaxWiener- 
Tree, and assume without loss of generahty that 2 < di < d2 < ■ ■ ■ < dn and that 
dn+i = ■ ■ ■ = dr = 1. (Since the case n = 1 is trivial, we assume from now on that 
n > 2.) It is straightforward to see that the backbone of the maximizing caterpillar will 
consist of n vertices vi, . . . ,Vn- In contrast to this, it is not straightforward to write 
down the sequence £i, . . . ,in of leaf-numbers: If one of the inner backbone vertices Vi 
with 2 < i < n — 1 gets degree di then it is adjacent to £i = di — 2 leaves, whereas if one 
of the two outermost backbone vertices Vi with i = 1 or i = n gets degree di then it is 
adjacent to £i = di — 1 leaves. 

Motivated by the discussion in the preceding paragraph, we introduce an (n + 2)- 
dimensional instance of the Wiener Max-QAP: the product matrix is built around the 
n numbers Oi = di — 1 for 1 < i < n and the two additional numbers On+i = = 1- 

The ID-distance matrix is built around the n points /3j = i — 1 for 2 < i < n + 1 and the 
two additional points /3i = 1 and /3n+2 = n. Furthermore, we impose the constraints 
that value an+i = 1 must be assigned to point /3i = 1, and that value a„+2 = 1 must 
be assigned to point (3n+2 = n. These additionally imposed constraints take care of the 
special treatment of the two outermost backbone vertices vi and f„. 

Now let us verify that the machinery of Sections [3] and S] still can be applied to this 
variant of the Wiener Max-QAP. First of all, the decomposition property works out 
exactly as before. Also Theorem 13.21 continues to hold, since adding 7r(l) = n + 1 and 
7r(n + 2) = n -|- 2 to a V-shaped permutation vr(2), . . . , 7r(n -|- 1) of the numbers 1, . . . , n 
always yields a V-shaped permutation. The dynamic program in Section [5] needs some 
small cosmetic changes that are caused by the additionally imposed constraints. 

• All states {k, m, L, R) with k < n must satisfy 2<m<n — k + 2 and L,R > 1. 

• All states {n + 1, m, L, R) must satisfy m = 1, L = 0, and R = 1. 

• All states (n -|- 2, m, L, R) must satisfy m = 1 and L = R = 0. 

These conditions ensure that the dynamic program assigns the values On+i and an+2 
during the last two stages to points /3i and Pn+2 exactly as desired. 

Furthermore, we note that in equations (jlSp and (I14p we always have |/3m+i —/3m I = 1 
and \/3m+k-i - Pm+k-2\ = 1 for /c < n, and |/3m+i - = and \(3m+k-i - Pm+k-2\ = 
\Pk — Pk-i\ for k > n. This implies that neither the recursive computations nor the 
values Z{k, m, L, R) do depend on the second coordinate m, which consequently may 
be dropped. (This should also be clear intuitively, since this second coordinate encodes 
the piece of the backbone to which the first k values ai,...,Qfc are assigned. All 
backbone pieces of length k are paths on k vertices and thus have the same combinatorial 
structure.) 

What about the time complexity? Exactly as in the proof of Theorem 14.11 the 
time complexity is proportional to the number of states {k,L,R). Since k can take 
0{r) possible values, and since L can take O(^Qj) = 0(r) possible values, the time 
complexity is O(r^). 
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Theorem 5.2 The problem of finding a tree that maximizes the Wiener index among 
all trees with a prescribed degree sequence can be solved in quadratic time O(r^), where 
r denotes the overall number of terms in the degree sequence. □ 



6 Conclusions 

We have introduced the Wiener Max-QAP, a special case of the quadratic assignment 
problem. We have provided a complete picture of the computational complexity of this 
special case: It is NP-hard in the ordinary sense, and it is solvable in pseudo-polynomial 
time. Our investigations also gave us a polynomial time algorithm for finding a tree 
that maximizes the Wiener index among all trees with a prescribed degree sequence, 
thereby settling a prominent open problem from chemical graph theory. 

One obvious open problem is to bring the quadratic time complexity O(r^) in The- 
orem [52] down to 0(r log r) or perhaps even down to linear time 0(r). 

Another open problem concerns the Wiener Min-QAP, where the goal is to minimize 
the objective value in ([H). It is an easy exercise to rewrite and to adapt the results of 
Sections [3] and S] to the minimization version: The minimization version always has an 
optimal solution that is pyramidal (which means that the permutation is first increasing 
up to some value i, and then decreasing from i onwards). And the minimization problem 
can be solved by dynamic programming in pseudo-polynomial time, within the same 
time complexity as that in Theorem 14. li The main gap in our knowledge concerns the 
complexity of the Wiener Min-QAP, and we pose the open problem of deciding whether 
it actually is NP-hard. 
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